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This  paper  investigated  methods  for  passive  estimation  of  the  bearing  to  a 
slowly  moving  acoustically  radiating  source.  The  mathematics  for  the  solution 
to  such  a problem  is  analogous  to  estimating  the  time  delay  (or  group  delay) 
between  two  time  series.  Since  the  estimation  of  time  delay  Is  Intimately  re- 
lated to  the  coherence  between  two  time  series,  a summary  of  the  properties  of 
coherence  is  presented. 
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The  term  coherence  has  several  different  meanings  and  indeed 

DEFINITIONS.  The  ONE  WE  USE  HERE  IS  THE  MAGNITUDE  SQUARED  OF  THE 
COEFFICIENT  OF  COHERENCY  DEFINED  BY  WEINER  IN  1930,  In  PARTICULAR, 
FOR  OUR  PURPOSES  HERE,  WE  DEFINE  THE  COHERENCE  BETWEEN  TWO 
STATIONARY  RANDOM  PROCESSES  A AND  B AS  THE  MAGNITUDE  SQUARED  OF 
THE  CROSS  POWER  SPECTRUM  DIVIDED  BY  THE  PRODUCT  OF  THE  TWO  AUTO 
POWER  SPECTRA.  The  COHERENCE  IS  A FUNCTION  OF  FREQUENCY  AND  HAS 
THE  USEFUL  PROPERTY  THAT  IT  LIES  BETWEEN  ZERO  AND  UNITY.  It  IS, 

IN  EFFECT,  A NORMALIZED  CROSS  SPECTRAL  DENSITY  THAT,  IN  SOME 
SENSE,  MEASURES  THE  EXTENT  TO  WHICH  TWO  RANDOM  PROCESSES  ARE 
SIMILAR.  For  EXAMPLE,  TWO  UNCORRELATED  RANDOM  PROCESSES  ARE  ALSO 
INCOHERENT;  THAT  IS,  THE  COHERENCE  IS  ZERO  BETWEEN  UNCORRELATED 

processes.  Further,  the  coherence  between  two  linearly  related 

PROCESSES  IS  UNITY. 


-Next  Slide  Please- 
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The  physical  problem  that  motivates  this  research  is  a desire 
to  passively  estimate  geographical  information  about  the  state 

OF  AN  ACOUSTIC  SOURCE.  In  THE  DEVELOPMENT  HERE  AN  ACOUSTIC  POINT 
SOURCE  RADIATES  SPHERICAL  WAVES,  RECEIVED,  FIRST,  AT  ONE  SENSOR 
AND  SOME  DELAYED  TIME  LATER,  AT  A SECOND  SENSOR.  THE  SOURCE  IS 
ASSUMED  STATIONARY  FOR  THE  OBSERVATION  PERIOD  AND  THE  SENSOR 
SEPARATION  IS  ASSUMED  KNOWN.  EACH  RECEIVED  WAVEFORM  IS  OBSERVED 
IN  THE  PRESENCE  OF  UNCORRELATED  NOISE.  THE  PROBLEM  WE  ADDRESS  IS 
HOW  TO  ESTIMATE  THE  TRAVEL  TIME  OF  THE  WAVEFRONT  OR  TIME  DELAY 
FROM  ONE  SENSOR  TO  THE  NEXT. 

The  importance  of  obtaining  a good  time  delay  estimate  is  that 

IT  CAN  BE  USED  TO  FIX  THE  SOURCE  LOCATION  ON  A HYPERBOLIC  LOCUS 
OF  POINTS.  OF  COURSE,  IF  WE  HAVE  THREE  SENSORS  WE  CAN  ESTIMATE 
TWO  TIME  DELAYS  AND  THE  INTERSECTING  HYPERBOLIC  CURVES  CAN  BE 
USED  TO  ESTIMATE  SOURCE  POSITION. 


-Next  Slide  Please- 
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GENERAL  CASE 


SPECIFIC  CASE 

X,(t)  - S(t)  ♦ n,(t) 

Xg(t)  - aS(t  + D)  + n2(t) 
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(Slide  3) 

In  the  general  case  we  can  model  the  acoustic  source  propaga- 
tion AND  NOISE  CORRUPTED  RECEPTION  AS  SHOWN  HERE.  In  PARTICULAR/ 
WE  TREAT  THE  PATH  FROM  THE  SOURCE  TO  EACH  RECEIVER  AS  A LINEAR 
TIME  INVARIANT  FILTER.  THE  RECEIVED  SIGNALS  X ONE  AND  X TWO 
CONSIST  OF  THE  FILTER  OUTPUTS  PLUS  NOISE. 

A SPECIAL  CASE  OF  THIS  MODEL  IS  SHOWN  ON  THE  BOTTOM  OF  THE 
SLIDE.  THE  FIRST  RECEIVED  WAVEFORM  CONSISTS  OF  SIGNAL  PLUS  NOISE. 
The  SECOND  RECEIVED  WAVEFORM  CONSISTS  OF  AN  ATTENUATED  AND  DELAYED 
SIGNAL  IN  THE  PRESENCE  OF  NOISE.  THE  MATHEMATICAL  PROBLEM  WE 
ADDRESS  IS:  HOW  TO  BEST  ESTIMATE  THE  TIME  DELAY  OR  EQUIVALENTLY 

SOURCE  BEARING.  FURTHER  WE  ARE  CONCERNED  WITH  THE  ROLE  OF  COHER- 
ENCE IN  THIS  PROCESS. 

For  analytic  purposes  we  treat  the  noise  as  stationary  and 

UNCORRELATED.  LATER  WE  MAKE  AN  IMPLICIT  ASSUMPTION  THAT  THE 
NOISE  IS  NORMAL  (GAUSSIAN). 


-Next  Slide  Please- 
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For  the  general  model  we  can  show  that  the  received  signal-to- 

NOISE  RATIO  IS  A FUNCTION  OF  ONLY  THE  COHERENCE  BETWEEN  THE  SOURCE 
AND  THE  RECEIVER.  INDEED,  RELATIVE  TO  THE  SENSOR  NOISE  POWER, 

THE  AMOUNT  OF  POWER  RECEIVED  FROM  THE  SOURCE  AFTER  IT  HAS  BEEN 
ATTENUATED  BY  ACOUSTIC  TRANSMISSION  THROUGH  THE  OCEAN  MEDIUM  IS 
DESCRIBED  BY  THE  SOURCE-TO-SENSOR  COHERENCE  DIVIDED  BY  ONE  MINUS 
SOURCE-TO-SENSOR  COHERENCE.  MOREOVER,  THE  COHERENCE  BETWEEN 
THE  TWO  RECEIVED  WAVEFORMS  CANfelQI  EXCEED  THE  COHERENCE  BETWEEN 
THE  SOURCE  AND  ANY  SENSOR;  THIS  IS  TRUE  WHEN  THE  OCEAN  MEDIUM 
IS  MODELED  AS  A LINEAR  TIME  INVARIANT  FILTER. 

More  specifically,  the  coherence  between  the  two  received 

WAVEFORMS  IS  EQUAL  TO  THE  PRODUCT  OF  THE  COHERENCES  BETWEEN  THE 
SOURCE  AND  EACH  OF  THE  RECEIVED  WAVEFORMS. 

-Next  Slide  Please- 
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Under  standard  assumptions,  namely,  that  the  two  received 

WAVEFORMS,  X ONE  AND  X TWO,  ARE  JOINTLY  STATIONARY,  NORMAL  (GAUSSIAN) 
RANDOM  PROCESSES  AND  THAT  THE  OBSERVATION  TIME  P IS  LARGE,  THE 
MAXIMUM  LIKELIHOOD  ESTIMATE  OF  TIME  DELAY  CAN  BE  DERIVED.  THE 
MAXIMUM  LIKELIHOOD, OR  ML, ESTIMATE  OF  TIME  DELAY  CAN  BE  INSTRUMENTED 
INONE  OF  TWO  WAYS.  SHOWN  HERE  IS  ONE  REALIZATION.  THE  FIRST 
RECEIVED  WAVEFORM  IS  FILTERED  BY  H ONE  TILDE,  AND  THE  SECOND 
RECEIVED  WAVEFORM  IS  FILTERED  BY  H TWO  TlLDE  AND  DELAYED.  THE 
FILTERS  MUST  HAVE  IDENTICAL  PHASE  RESPONSES.  THE  FILTER  OUTPUTS 
ARE  SUMMED,  SQUARED,  AND  INTEGRATED  AS  SHOWN.  THE  HYPOTHESIZED 
VARIABLE  DELAY  THAT  MAXIMIZES  THIS  SYSTEM  OUTPUT  IS  THE  MAXIMUM 
LIKELIHOOD  ESTIMATE  OF  TIME  DELAY.  THE  SPECIFIC  FILTER  CHARAC- 
TERISTICS DEPEND  ON  THE  SIGNAL  AND  NOISE  SPECTRA,  WHICH  MUST  BE 
KNOWN  OR  ESTIMATED. 

-Next  Slide  Please- 
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Another  realization  for  the  maximum  likelihood  estimate  of 

TIME  DELAY  IS  A SPECIAL  CASE  OF  THE  GENERALIZED  CROSSCORRELATION 
PROCESSOR.  In  THIS  PROCESSOR  THE  FIRST  RECEIVED  WAVEFORM  IS  FIL- 
TERED BY  H ONE  AND  THE  SECOND  RECEIVED  WAVEFORM  IS  FILTERED  BY 
H TWO,  DELAYED,  MULTIPLIED,  AND  INTEGRATED  AS  SHOWN  ON  THE  DIAGRAM 
AT  THE  TOP  OF  THE  SLIDE.  By  PROPER  CHOICE  OF  THE  GENERAL  WEIGHTING 
FUNCTION,  W,  WHICH  IS  THE  PRODUCT  OF  H ONE  AND  H TWO  CONJUGATE, 

WE  CAN  ACHIEVE  THE  MAXIMUM  LIKELIHOOD  ESTIMATE  FOR  TIME  DELAY. 

However,  we  can  also  achieve  a generalized  crosscorrelation 

FUNCTION,  R OF  TAU,  BY  MULTIPLYING  THE  CROSS  SPECTRUM,  G, 

BETWEEN  THE  TWO  RECEIVED  WAVEFORMS,  BY  THE  GENERAL  WEIGHTING,  W, 

AND  COMPUTING  THE  FOURIER  TRANSFORM  AS  INDICATED  ON  THE  BOTTOM  OF 
THE  SLIDE. 

-Next  Slide  Please- 
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If  we  estimated  the  generalized  crosscorrelation  function 

FOR  SIX  DIFFERENT  TRIALS,  THE  PEAK  OF  THE  FUNCTION  MIGHT  VARY 
AS  A FUNCTION  OF  TRIAL.  We  HAVE  ACTUALLY  IMPLEMENTED  THE  TECHNI- 
QUE FOR  SEVERAL  EXAMPLE  CASES  ON  THE  UnIVAC  1108.  BASED  ON  OUR 
EXPERIMENTAL  RESULTS,  WE  SPECULATE  THAT  A TYPICAL  GENERALIZED 
CROSSCORRELATION  FUNCTION  MIGHT  PEAK,  AS  INDICATED  IN  THE  HYPO- 
THETICAL TRIALS  SKETCHED  HERE.  In  PARTICULAR,  THE  ABSCISSA  VALUE 
OF  THE  PEAK  LOCATION,  THAT  IS,  THE  ESTIMATE  OF  TIME  DELAY,  HAS  A 
CERTAIN  AMOUNT  OF  VARIATION.  NOTICE  ALSO  IN  TRIAL  NUMBER  5,  NEXT 
TO  THE  BOTTOM  PLOT,  THAT  A NUMBER  OF  AMBIGUOUS  PEAKS  CAN  ARISE  IN 
ADDITION  TO  THE  LOCAL  VARIATION  OF  THE  TIME  DELAY  ESTIMATE.  THE 
AMBIGUITY  PROBLEM  IS  NOT  TREATED  IN  THIS  WORK.  THE  PROBLEM  OF 
COMPUTING  THE  VARIANCE  OF  THE  TIME  DELAY  ESTIMATE  IS  A DIFFICULT  • 
ONE  IN  WHICH  ONE  IS  PUZZLED  HOW  TO  PROCEED.  HOWEVER,  IF  WE  COULD 
COUNT  THE  NUMBER  OF  PEAKS  THAT  OCCURRED  AT  EACH  OF  SEVERAL 
ABSCISSA  VALUES,  THEN  WE  COULD  PLOT  A FREQUENCY  DISTRIBUTION  OF 
THE  PEAK  LOCATION.  FROM  THIS  DISTRIBUTION  WE  CAN  OBTAIN  THE 
VARIANCE  OF  THE  TIME  DELAY  ESTIMATE. 

-Next  Slide  Please- 
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The  variance  of  the  time  delay  estimate  is  a complicated 

FUNCTION  OF  SEVERAL  PARAMETERS.  It  DEPENDS  ON  THE  LENGTH  OF 
THE  OBSERVATION  TIME/  ?,  THE  GENERAL  WEIGHTING  FUNCTION/  W/  THE 
AUTO-SPECTRAL  DENSITIES  OF  THE  TWO  RECEIVED  WAVEFORMS/  AND  THE 
MAGNITUDE  CROSS  SPECTRUM  BETWEEN  THE  TWO  RECEIVED  WAVEFORMS. 

It  ALSO  DEPENDS  ON  THE  COHERENCE/  C/  DEFINED  earlier  as  the 
MAGNITUDE  SQUARED  CROSS  SPECTRUM  DIVIDED  BY  THE  PRODUCT  OF  THE 
TWO  AUTO-SPECTRAL  DENSITIES.  RECALL  THE  COHERENCE  IS  GREATER 
THAN  OR  EQUAL  TO  ZERO  AND  IS  LESS  THAN  OR  EQUAL  TO  UNITY.  WHEN 
THE  OBSERVATION  TIME  IS  LARGE  OR  THE  COHERENCE  IS  NEAR  UNITY/  THE 
VARIANCE  IS  GENERALLY  QUITE  LOW  AND  YOU  CAN  DO  WELL  IN  SPITE  OF 
THE  WEIGHTING  SELECTED.  Of  COURSE/  AN  IMPORTANT  ROLE  TO  BE 
PLAYED  BY  THE  EXPRESSION  HERE  IS  TO  EVALUATE  HOW  DIFFERENT  PRO- 
CESSORS COMPARE  WITH  ONE  ANOTHER.  ANOTHER  IMPORTANT  USE  OF 
THIS  EXPRESSION  IS  IF  ONE  KNOWS  THEORETICALLY  THE  BEST  WEIGHTING 
FUNCTION  TO  APPLY/  BUT  APPLIES  AN  INCORRECT  OR  SUBOPTIMUM  WEIGHT- 
ING/ THEN  THE  VARIANCE  OF  THE  SUBOPTIMUM  DELAY  ESTIMATOR  CAN  BE 
EVALUATED. 


-Next  Slide  Please- 


17 


TO  5507 


1 


MINIMUM  VARIANCE 


2P  f df  (2trf) 


2_CliW_ 

1-C12(f) 


-1 


FOR 


WU1  (f)  = 


C12(f) 


ML 


g,2^)I  Fi  “ c 


.*«] 


Slide  9 


18 


(Slide  9) 


TD  5507 


The  minimum  variance  for  any  time  delay  estimation  scheme 

CAN  BE  OBTAINED  FROM  THE  CRAMER  RaO  LOWER  BOUND.  As  SHOWN  HERE/ 

IT  IS  A FUNCTION  OF  ONLY  TWO  PARAMETERS:  THE  OBSERVATION  TIME  P 

AND  THE  COHERENCE  BETWEEN  THE  TWO  RECEIVED  WAVEFORMS.  As  P IS 
INCREASED  THE  VARIANCE  DROPS;  FURTHER/  AS  THE  COHERENCE  C TENDS 
TOWARD  UNITY  THE  TERM  C OVER  ONE  MINUS  C SQUARED  TENDS  TOWARDS 
INFINITY.  THUS/  AS  THE  COHERENCE  OR  C TENDS  TOWARDS  UNITY/  THE 
VARIANCE  TENDS  TOWARDS  ZERO.  HOWEVER/  THE  COHERENCE  IS  NOT  UNDER 
OUR  CONTROL.  The  FACTORS  WHICH  WE  CAN  CONTROL  ARE  THE  OBSERVATION 
TIME  P AND  THE  WEIGHTING.  THE  MINIMUM  VARIANCE  IS  ACHIEVED  FOR 
THE  MAXIMUM  LIKELIHOOD  WEIGHTING  FUNCTION  GIVEN  BY  C OVER  ONE 
MINUS  C TIMES  THE  MAGNITUDE  CROSS  SPECTRUM. 

-Next  Slide  Please- 
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The  maximum  likelihood  weighting  multiplies  the  estimated 

CROSS  SPECTRUM  TO  YIELD  A SINGLE  FUNCTION  TO  BE  FOURIER  TRANS- 
FORMED. In  general,  when  the  true  values  OF  COHERENCE  AND 
MAGNITUDE  CROSS  SPECTRUM  ARE  UNKNOWN,  THEY  MUST  BE  ESTIMATED. 

Estimates  are  indicated  by  hats. When  spectral  analysis  is  used 

TO  YIELD  ESTIMATES  IN  PLACE  OF  THE  TRUE  QUANTITIESj  THE  FUNCTION 

to  be  Fourier  transformed  is  indicated  on  this  slide-  The  cross 

SPECTRUM  OVER  THE  MAGNITUDE  CROSS  SPECTRUM  CAN  BE  THOUGHT  OF  AS 
E TO  THE  MINUS  J PHASE.  In  PARTICULAR,  NOTE  THAT  THE  WEIGHTING 
EMPHASIZES  THE  PHASE  OF  THE  ESTIMATED  CROSS  SPECTRUM  IN  THOSE 
FREQUENCY  BANDS  WHERE  THE  COHERENCE  IS  HIGH.  ONE  WOULD  EXPECT 
THE  ESTIMATED  PHASE  OF  THE  CROSS  SPECTRUM  TO  PLAY  AN  IMPORTANT 
ROLE  IN  TIME  DELAY  ESTIMATION,  SINCE  THE  SLOPE  OF  THE  PHASE  IS 
A MEASURE  OF  THE  TIME  DELAY.  We  CAN  SEE  THIS  BY  NOTING  THAT  THE 
PHASE  SLOPE  IS  MEASURED  IN  RADIANS,  DIVIDED  BY  RADIANS  PER  SECOND, 
OR  SECONDS.  OF  COURSE,  THE  PHASE  ESTIMATES  WILL  BE  NOISY  IN  THOSE 
FREQUENCY  BANDS  WHERE  THE  COHERENCE  IS  LOW  SO  WE  WILL  EMPHASIZE 
THE  PHASE  IN  THOSE  BANDS  WHERE  THE  COHERENCE  IS  HIGH. 

-Next  Slide  Please- 
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SUMMARY 

• ACOUSTIC  SOURCE 

• TIME  DELAY  MODEL 

• DERIVED  ML  TIME  DELAY  ESTIMATE 

• DERIVED  CRAMER  RAO  LOWER  ROUND 

• DERIVED  THE  VARIANCE  FOR  ANY  GCC 

• SHOWN  ML  ESTIMATE  IS  MINIMUM  VAR 

• IMPLEMENTED  RESULTS 

• APPLICATIONS  TO  ESTIMATING  SOURCE  POSITION 
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In  summary,  the  physical  problem  motivating  this  research  is 

A DESIRE  TO  ESTIMATE  POSITIONAL  INFORMATION  ABOUT  AN  ACOUSTIC 
SOURCE.  WE  HAVE  PROPOSED  A TIME  DELAY  MODEL  AND  DERIVED  THE 
MAXIMUM  LIKELIHOOD  ESTIMATE  FOR  TIME  DELAY.  ADDITIONALLY  WE  HAVE 
DERIVED  THE  CRAMER  RAO  LOWER  BOUND  ON  THE  VARIANCE  OF  THE  TIME 
DELAY  ESTIMATE.  SUBSEQUENTLY  WE  HAVE  DERIVED  AN  EXPRESSION  FOR 
THE  VARIANCE  OF  THE  TIME  DELAY  ESTIMATE  FOR  ANY  GENERALIZED  CROSS- 
CORRELATION PROCESSOR.  WE  HAVE  SHOWN  THAT  THE  MAXIMUM  LIKELIHOOD 
ESTIMATE  OF  TIME  DELAY  ACHIEVES  THE  CRAMER  RAO  LOWER  BOUND  AND  IS 
THEREFORE  MINIMUM  VARIANCE;  AS  SUCH  THE  PROPOSED  TECHNIQUE  IS  THE 
BEST  PROCESSING  THAT  CAN  BE  DONE  TO  ESTIMATE  TIME  DELAY  OR.  EQUIV- 
ALENTLY, TO  ESTIMATE  THE  HYPERBOLIC  LOCUS  OF  POINTS  ON  WHICH  THE 
ACOUSTIC  SOURCE  IS  LOCATED.  THERE  IS  NO  BETTER  TECHNIQUE.  We 
HAVE  IMPLEMENTED  THE  RESULTS  IN  AN  APPROXIMATE  METHOD  BY  SUB- 
STITUTING ESTIMATED  MAXIMUM  LIKELIHOOD  WEIGHTING  IN  PLACE  OF 
TRUE  WEIGHTING  AND  FOUND  THAT  THE  TECHNIQUE  WORKS  ON  A LARGE 
SCALE  DIGITAL  COMPUTER.  Of  COURSE,  THE  ABILITY  TO  LOCATE  A 
SOURCE  ON  A HYPERBOLIC  LOCUS  OF  POINTS  SUGGESTS  THAT,  WITH  THREE 
SENSORS,  INTERSECTING  HYPERBOLIC  CURVES  CAN  BE  USED  TO  ESTIMATE 


SOURCE  POSITION. 
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THE  ROLE  OF  COHERENCE  IN  TIME  DELAY  ESTIMATION 

G.  Clifford  Carter 

Naval  Underwater  Systems  Center 
New  London,  Connecticut  06320  U.S.A. 


ABSTRACT.  This  paper  investigates  methods  for  passive  estimation 
of  the  bearing  to  a slowly  moving  acoustically  radiating  source. 
The  mathematics  for  the  solution  to  such  a problem  Is  analogous 
to  estimating  the  time  delay  (or  group  delay)  between  two  time 
series.  Since  the  estimation  of  time  delay  is  intimately  related 
to  the  coherence  between  two  time  series,  a summary  of  the  pro- 
perties of  coherence  is  presented. 

The  maximum  likelihood  (ML)  estimate  of  time  delay  (under 
jointly  stationary  Gaussian  assumptions)  is  presented.  The 
explicit  dependence  of  time  delay  estimates  on  coherence  is  evi- 
dent in  the  estimator  realization  in  which  the  two  time  series 
are  prefiltered  (to  accentuate  frequency  bands  according  to  the 
strength  of  the  coherence)  and  subsequently  crosscorrelated.  The 
hypothesized  delay  at  which  the  generalized  crosscorrelation  (GCC) 
function  peaks  is  the  time  delay  estimate.  The  variance  of  the 
time  delay  estimate  is  presented  and  discussed. 


INTRODUCTION.  An  acoustic  source  whose  signal,  s(t),  is  trans- 
mitted through  the  ocean  medium  and  received  in  the  presence  of 
additive  noise  can  be  characterized  by 


x^t)  = s i ( t ) 


+ n.(t) 


1 ■ 1,2 


(1) 


For  the  main  purposes  of  this  paper  s.(t)=s(t),  S2(t)*as(t+D)  and 
we  desire  to  present  an  ML  estimator  *for  the  time  delay  D.  Tne 
delay  parameter  can  be  used,  in  a nondispersi ve  medium  with  known 
speed  of  transmission,  to  estimate  the  bearing  to  an  acoustic 
source  (relative  to  the  sensor  baseline)  or,  more  generally,  to 
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estimate  a hyperbolic  "line"  of  position.  Since  the  final  result 
depends  heavily  on  the  coherence  between  xj  and  X2,  we  precede 
the  development  with  a concise  review  of  the  properties  of  the 
coherence  function  and  of  results  that  bear  directly  on  the  esti- 
mation of  time  delay. 

THEORY  OF  COHERENCE.  For  any  two  jointly  stationary  random  pro- 
cesses xi  and  X2,  the  coefficient  of  coherency  or  the  complex 
coherence  has  been  defined  by  Weiner  (1930)  as  the  ratio 


where  Gv  (f)  is  the  cross  power  spectral  density  function  be- 
xlx2 

tween  xj  and  X2,  and  G (f),  1*1,2  are  the  auto  power  spectral 

xixi 

density  functions  at  frequency,  f. 

The  magnitude-squared  coherence  (MSC)  or  simply  the  coherence 
is  defined  by  (see,  for  example.  Carter,  Knapp  and  Nuttall  (1973)) 


A useful  property  of  the  MSC  is 


0 < C 


(f)  £ 1 


provided  the  autospectra  are  positive  (in  particular  non  zero). 

In  order  to  attach  some  physical  significance  to  what  the 
coherence  measures,  consider  that  the  ocean  medium  operators  Mi 
and  M2  are  linear  time-invariant  filters.  Thus  Sj(t)  and  S2(t) 
in  equation  (1)  are  the  respective  outputs  of  filters  Mj(f)  and 
M2(f)  when  excited  by  source  s(t).  When  the  noise,  nj(t),  is  un- 
correlated with  the  signal,  s(t),  at  the  i-th  sensor,  the  ratio 
of  the  received  signal  power  at  the  output  of  the  ocean  channel 
to  the  corruptive  noise  power  depends  on  the  coherence  between 
the  source  and  the  sensor.  Specifically,  from  Carter,  Knapp, 
and  Nuttall  (1973) 
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6s  s (f)  Csx  (f) 
S1S1  _ SX1 


Vi(f)  1'S(f) 


1 - 1,2 


(3) 


That  Is,  the  received  signal-to-noise  ratio  (SNR)  at  the  i-th sen- 
sor depends  on  the  coherence  between  the  source  and  the  received 
waveform.  This  result  has  been  expressed  by  Carter  and  Knapp 
(1976)  more  compactly  as 


Cx  x <f> 
XjX2 


C~  <f>  Csx.'fl 


Sx 


(4) 


These  results  apply  only  to  the  case  where  the  medium  can  be 
accurately  modeled  by  linear  time-invariant  filters  corrupted  by 
uncorrelated  additive  noise. 


RESULTS.  For  the  purpose  of  obtaining  an  ML  estimate  of  delay, 
certain  assumptions  are  required.  In  particular,  for  a signal 
emanating  from  a nearfield  source  and  monitored  in  the  presence  of 
noise  at  two  spatially  separated  sensors  we  require  in  equation  (1) 
that  Sj(t)  = s(t)  and  S2 ( t ) = as(t+D).  Further,  we  require  that  a 
is  real  and  s(t),  nj(t),  and  n^(t)  are  real,  jointly  stationary, 
Gaussian  random  processes.  Source  s(t)  and  noises,  (t ) andn2(t) 
are  assumed  to  be  mutually  uncorrelated. 


An  estimated  value  of  D is  the  hypothesized  value  t that  maxi- 
mizes the  generalized  crosscorrelation  (GCC)  function  defined  by 


R(t) 


■L 


xlx2 


(f)W(f)e 


j2irfx 


df. 


(5) 


For  xi(t)  and  xg( t)  real,  the  ML  estimator  requires  a particular 
hting. 


weighti 

W(f)  - Hj(f)H2*(f)  = 


xlx2 


(f) 


lGx,»2(,)l  [’  - cv2(,)] 


(6) 


A complete  derivation  is  given  by  Carter  (1976). 


Note  from  equation  (6)  that  for  the  ML  estimate  of  delay  that 
W(f)  is  real.  The  ML  estimator  is  virtually  equivalent  to  one 
proposed  by  Hannan  and  Thomson  (1973).  The  ML  estimator  can  be 


achieved  by  shaping  Xj(t)  with  filter  H.(f)  and  x2(t)  with  filter 
Hgff)  crosscorrelating  the  filter  outputs,  and  observing  what  hypo- 
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theslzed  value  of  delay  achieves  a maximum. 

The  estimator  can  also  be  achieved  by  other  methods.  For 
example,  Hahn  (1975),  Carter  and  Knapp  (1976)  and  Carter  (1976) 
present  a method  of  filtering  and  sunning  the  outputs,  squaring 
and  averaging  in  order  to  estimate  the  delay  D.  The  processor 
could  also  be  realized  as  a number  of  "best"  estimates  of  D for  a 
variety  of  frequencies.  The  ML  estimate  Is  then  achieved  by  per- 
forming a weighted  average  across  frequency.  For  example.  Clay, 
Hinich  and  Shaman  (1973)  develop  ML  estimates  of  bearing  (analo- 
gous to  delay)  for  each  of  a number  of  different  frequencies.  To 
obtain  a single  estimate  of  source  bearing,  these  individual  esti- 
mates should  then  be  combined  with  weighting  dependent  upon  the 
■particular  underlying  signal  and  noise  characteristics. 

The  role  of  coherence  in  the  weighting  used  for  ML  estimation 
of  D is  specified  in  equation  (6).  Note  that  those  values  of  co- 
herence near  unity  are  most  important;  conversely,  in  those  fre- 
quency bands  where  there  is  no  source  signal  power  (hence,  where 
the  received  waveforms  are  incoherent),  the  delay  estimate,  as 
would  be  expected,  receives  no  weight.  The  ML  estimator  is  actu- 
ally a function  of  more  fundamental  spectral  measurements  than 
those  specified  in  equation  (6).  However,  expressing  the  pro- 
cessor In  more  fundamental  but  unnormalized  quantities  can  make 
Interpretation  more  difficult,  though  equally  correct. 

The  ML  weighting  agrees  with  MacDonald  and  Schulthelss 
(1969),  and  Hahn  (1975)  under  specific  conditions  (including  when 
there  are  two  sensors  and  no  attenuation). 

VARIANCE  OF  GENERAL  TIME  DELAY  ESTIMATORS.  The  variance  of  the 
time  delay  estimate  in  the  neighborhood  of  the  true  delay  for 
general  weighting  function  W(f)  is  given  by 


Varjp-p)=  f?{fi  [‘  - Cl?(f)]df  (7) 


‘ /• 00 
/Wr 

•/-CD 


lGx  x 
Xlx2 


W(f)df  2 


where  P 
(6)  and 


is  the  observation  period  (in  seconds). 
(7),  the  variance  of  the  ML  processor  is 


VarML  [d-d|= 

’ r (2»f)2  c12(f) 

2P  / l£  df 

L J 

0 1 - C12<f>  J 

From  equations 


1 


(8) 
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The  ML  processor  achieves  the  Cramer-Rao  lower  bound  (see  Carter 
(1976)).  Therefore,  the  ML  processor  achieves  a variance  less 
than  or  equal  to  that  provided  by  other  correlation  processors. 

These  results  for  variance  can  be  related  to  MacDonald  and 
Schultheiss  ( 1969)  as  follows.  Define  the  bearing  to  an  acoustic 
source,  as  in  Nuttall,  Carter  and  Montavon  (1974) 

4 * arc  cos  ^ (9) 

where  £ is  the  speed  of  sound  in  the  nondispersive  medium  and  d 
is  the  sensor  separation.  Consider  the  case  where  the  estimated 
D equals  the  true  delay  plus  a perturbation.  By  a Taylor  series 
expansion, it  follows  for  the  bearing  error  defined  by  the  differ- 
ence between  the  true  bearing  and  the  estimated  bearing  that  the 
standard  deviation  of  the  bearing  error  is  given  by  (Carter 
(1976)): 


[v«r  (♦-♦)]  S [var  ( D - 0 )] 


The  term  d sin  $ can  be  viewed  as  the  effective  array  length 
(sensor  separation)  physically  steered  at  the  source. 

The  combining  of  equations  (8)  and  (10) 
suggests  that,  in  order  to  reduce  the  variance  of  the  bearing 
estimate,  the  observation  period  and  the  sensor  separation  should 
be  made  as  large  as  possible.  This  agrees  with  one's  intuition 
and  the  results  of  MacDonald  and  Schultheiss  (1969).  Further, 
the  fact  that  equation  (10)  depends  on  the  effective  array  length 
physically  steeled  toward  the  source  suggests  the  desirability 
of  sensor  mobility  to  maximize  sin  $ when  d is  limited. 
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